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Abstract 


In this study, a model was developed for predicting the seeding rate of corn seeds of a typical row- 
crop planter equipped with a multi-slot feeding device. To this, nine multi-slot rotors (with 4, 5 and 6 
slots in three angles of mouth including 23°, 25° and 27°) were designed and manufactured. Tests 
were carried out at four levels of angular velocity of 40, 52, 62 and 78 rpm on grease belt moving at 
constant speed of 3.5 km h’'. Tests were completed in three replications. Discharge flow rate was 
measured and recorded for each treatment. The data were used to develop a model which can be used 
for predicting the seeding rate under various numbers of slot, mouth angle and rotor angular velocity. 
According to the results, angle mouth of slots, number of slots, angular velocity and the dual 
interaction between them showed increasing effects on weight flow rate of seeds (P-value<0.01). In 
the next step, raw data were used to develop the two desired models: based on the dimensional 
analysis technique and response surface methodology (RSM). The models outputs were compared to 
experimental data. The standard error of estimate for flow rate for dimensional analysis and response 
surface methodology (RSM) were 68.13 mm’ s' and 475.59 mm’ s", respectively. The dimensional 
analysis model was closer to experimental data rather than the RSM method. Thus, to predict the 
volume flow rate of seed, the dimensional analysis model is recommended. 


Keywords: Corn, Dimensional analysis, Response surface methodology, Row planting 


Notation 
Vertical intersept C_ Solids volumetric Flowrate (mm* s') Q 
Effective diameter of seed (mm) d Number of rotor R 
Effective diameter of Rotor(mm) D, Rotor thikness (mm) t 
Casing inside diameter(mm) D_ The mean width of slots (mm) W 
Acceleration of gravity (=9806) g Solid mass flow rate (g min’) W 
Depth of the slot (mm) h_ Angle of repose (degree) 8 
Slope of logarithmic line K Angle mouth of slot (degree) B 
Mean seed longitude (mm) 1 Density of solid particles (g min’) Ds 
Rotating speed of rotor (rps) n Loose bulk density (g min”) y 
Slots number N_ Function of.... Q 
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Introduction 


Some of the seeds planted per unit area 
cannot survive through germination and 
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emergence stages. Hence, the seeds number 
must be greater than the desired plant stand 
(Raheman and Singh, 2003; Srivastava, 2006). 
Therefore, a planter must meter seeds at proper 
rate and also must control the horizontal 
placement of seeds in the desired pattern. 

The seed metering mechanism is one of the 
main components of a planter that affects 
directly on yield efficiency (Ani et al., 2016). 
Seed metering devices are categorized as 
either precision or dune flow metering based 
on planting pattern result (Murray ef al., 
2006). Single-seed metering devices have been 
developed to provide precision for seed 
spacing planters which can provide exact and 
unique placement of single seeds at equal 
intervals within each rows. However, in this 
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technique, less attention has been given to the 
physical properties of seeds (Khan, 2008). 

Precision meters attempt to meter single 
seeds. While there is a large range of precision 
metering devices, most of them can be broadly 
classified as ‘plate’, ‘belt’, ‘disc’, ‘drum’ or 
‘finger’ types. The classification largely 
depends on the design and/or shape of the 
principle moving element that enables seed 
simulation (i.e. the selection of single seeds 
from the seed lot) (Gray and MacIntyre, 2012). 
Plate metering devices consist of a driving 
disk with indents on its periphery, and their 
efficiency depends on the seeds size (Ani et 
al., 2016). They can be sub-categorized as 
horizontal and vertical plate types (Murray et 
al., 2006). In this context, the multi-slot rotors 
are classified among the vertical sheet metal 
metering devices. 

Among the agricultural crops, corn has 
great importance as food, livestock and poultry 
feed, industrial and medicinal consumptions 
and as a rich source of carbohydrates, protein, 
iron, and vitamin B (Al-Mallahi and Kataoka, 
2013; Ani et al., 2016). It is necessary to 
cultivate corn seeds in a row. Therefore, the 
distance between consecutive seeds should be 
kept constant to minimize damage to them. 
This distance depends on the correct forward 
speed of the planter, rotational speed of the 
metering device and the space of the metering 
device cells (Anantachar et al., 2011). 

Rebati and Zareian (2003) fabricated a 
metering device for hill-planting with parallel 
slots on its roller surface which was used for 
hill-planting of roof rice in _ laboratory 
conditions. In order to evaluate the 
performance of the unit, forward velocity and 
diameter of the seed metering device were 
considered. According to their findings, the 
seeding speed of 3 km h” was considered as 
the optimum speed for minimum seeds’ 
damage. In another investigation, Jafari (1991) 
considered main characteristics of a fluted- 
roller metering device of a grain drill and 
developed a model equation for predicting 
discharge volume flow rate. 

A new methodology to estimate the mass of 
grain seeds, which flow in the shape of 


clumps, was suggested in a study. The 
methodology used an off-the-shelf digital fiber 
sensor to detect the behavior of the clumps and 
multiple linear regression modeling to estimate 
the mass by the parameters detected by the 
sensor which were the length and the density 
of the clumps. An indoor apparatus was used 
for modeling which resembled the sowing 
process using the grain drill. A fluted roller 
was installed in the apparatus to regulate the 
flow of seeds. Results showed that the digital 
fiber sensor could be used for estimating mass 
flow of seeds at variable sowing rates within 
the speed limits of the grain drill (Al-Mallahi 
and Kataoka, 2013). 

In the absence of any _ characteristic 
relationship for predicting volume flow rate of 
a multi-slot rotors in row-crop planting, the 
current investigation was aimed to fabricate a 
test rig and establish a model for estimating 
the volume flow rate of corn seeds. 


Materials and Methods 


Laboratory tests were conducted to measure 
the discharge flow rate of feeding device for a 
number of treatments. Treatments were three 
rotors having different slots (4, 5 and 6), each 
rotor having three angle of mouth (23°, 25° 
and 27°) at four angular velocity (40, 52, 63 
and 78 rpm). Tests were replicated three times. 
Development of the required test-rig 

A test rig was assembled to measure the 
discharge flow rate for each settings (Figure 
1), the rig consists following parts: 

1- Metering device: In the design of the 
metering device, the shape and number of slots 
are important factors. The dimensions of each 
slot are associated with the seed size and the 
number of seeds cells. Three fundamental 
design variables including the effective 
diameter of the metering devices (D-), the slot 
depth (h) and the angled mouth of slots (8), as 
shown in Figure 2 (Rebati and Zareian, 2003). 
Considering the shape and size of the corn 
seeds, three angles of 23°, 25° and 27° were 
calculated and defined as the angles mouth of 
slots. Angles have direct effect on the depth of 
the slots and the effective diameter of the 
metering devices. In order to provide different 
distance between consecutive seeds on the 
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ground, the fluted rotors with 4, 5 and 6 slots Dassault Systemes,  Vélizy-Villacoublay, 
were used. Therefore, the nine rotors were France, 2016) and were fabricated from 
designed in SolidWorks software (SP 3.0, Polytetrafluoroethylene material (Figure 3). 


Fig.2. Variables related to design of the metering devices 


Fig.3. Metering devices used in this study 


2- Metering device Case: The metering material with dimensions of 50 mmx100 
device case was fabricated from polyethylene mmx75 mm. This case has a 50 mmx25 mm 
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intake gate and 50 mmx13 mm outlet opening, 
and a ball bearing (B17-99) was used to 
accommodate the shaft of the metering device 
in the center of the case. 

3- Power transmission system: In order to 
facilitate the movement of the metering device 
and control the spacing of seeds along the row 
of the grease belt, a chain, and sprocket system 
was used to provide the desired rotational 
speeds of 40, 52, 63 and 78 rpm. 

Experimental Design 

The volumetric flow rate of seed discharge 
from a planter feeding device depends on the 
following groups of parameters (Jafari, 1991): 

1. Metering device parameters, which 
include the effective diameter of the metering 
devices, the number of slots, the angled mouth 
of slots, the depth of slots and the rotational 
speed of the metering device. 

2. Seed parameters including mean length 
and the effective diameter of seeds, the density 
of solid particle, surface characteristics, and 
seed shape factor. 

3. Flow parameters including solids flow 
rate, gravitational acceleration, and friction 
between solids and wall. 

The dimensional analysis and response 
surface methodology (RSM) approaches were 
considered to express the relationship between 
the selected parameters. In summary, the bulk 
volume flow rate of solids (depended variable) 
was considered as a function of nine 
independent variables as follow: 

Q= (1,t,g,h,N, De, B,1, d) (1) 

Thus, based on z- Buckingham theorem 
(Murphy, 1950), the selected model has the 
formula for the linear form: 
nt? -1 
=C [o: (en te) X Pz (en e*) X @3(th-*) (2) 
x @4(NtDe~*) x @5(B) x @s(Del™) x @7(Ded)] 

The Response Surface Methodology (RSM) 
consists of a group of mathematical and 
Statistical techniques used in the development 
of an adequate functional relationship between 
a response of dependent and some independent 
variables denoted (Maged ef al., 2018; Khuri, 
2017; Salar and Karparvarfard, 2017). The 
miss index, multiple index, quality of feed 
index and precision in spacing were 


considered as dependent variables as details 
mentioned in Balanian ef al. (2018). The angle 
mouth of slots (§), number of slots (N) and the 
rotational speed of the rotor (n) are used as 
independent and changeable variables, and the 
volume flow rate of solids (Q) was considered 
as the dependent variable. The RSM model 
was recommended as below: 
Q=C,xB+C,xN+C3;xn+C,xBxN 
+C,xBxn+C,xNxn (3) 
+C, x B?+ Cg x N*+Cy 
xn? +C 
Where C;, Co, ..., Cg are constants. 
Considering the proposed models, the 
parameters in the study were selected. The 
average particle diameter and the length of 
scorn seeds were calculated according to a 
previously published method (Mohsenin, 
1970). Besides, the loose bulk density was 
calculated by using the equation (4): 
M 


Ve (4) 


Where M is the weight of the sample and V 
is the volume of the particle samples. The 
loose bulk density as obtained to be 0.81 g.L'! 
which was used to calculate the volume flow 
rate of particles. The thickness and outside 
diameter of each rotor were considered to be 
constant and equal to 15 mm and 60 mm, 
respectively. 

The number of slots, the depth of slots, the 
angled mouth of slots and the effective 
diameter of the rotors could be varied. By 
changing the values of the slot’s mouth angle, 
the quantities of h and D,. were defined by 
using SolidWorks® software. The 
specification and the weight flow rate of solid 
for all rotors are given in Table 1. Based on 
the size of the rotor case chamber, the rotor 
thickness was considered as 18 mm. At 
different times, the average values for weight 
flow rate (W) were obtained in three 
replications for each run. The forward speed of 
the grease belt was considered to be constant 
and adjusted as 3.5 km h''. Particles damage 
within the range for this test run were found 
approximately 5.6%. 
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Table 1- The specification and weight flow rate of Corn seed for each rotor used 


Rotor Number N f (degree) h(mm) D,(mm) n(rpm) W(g min’) 
40 88 
R, 4 23 7.72 52.22 ae L089 
63 122.83 
78 124.53 
40 93.55 
2 116.17 
Ro 5 23 3 52.22 2 bees 
78 136.40 
40 138.37 
52 180.94 
R; 6 23 7.72 52.22 me ee 
78 211.32 
40 97.43 
Ry 4 25 872 51.22 oe ere 
78 143.85 
40 125.46 
52 163.41 
Rs 5 25 8.72 51.22 63 181.60 
78 202.96 
40 157.43 
52 187 
Re 6 25 8.72 S122 Pe ios ee 
78 253.11 
40 120.85 
LY 149.82 
R 4 om 9.72 50.19 
: 63 180.98 
78 194.97 
40 148.79 
52 184.09 
Re 5 P| 9.72 50.19 ee Dae 
78 233.89 
40 191.81 
Ry 6 27 9.72 50.19 a mar 
78 328.91 


The regression equation allows us to make 
predictions, while no indication of the 
accuracy of the prediction was detected. Thus, 
the standard error of the estimate as the 
measure of the accuracy of prediction was 
calculated by the following equation: 

The standard error of estimate 


iP (measured value — predicted value)? 


(5) 


the number of pairs of scores 


Results and Discussion 


The Experimental Results 
By increasing the rotational speed of rotors, 
the weight of solids would be increased. This 


increase showed a descending slope which 
may be correlated with further increasing the 
centrifugal force and decreasing the seed 
entrance rate to the intake gate (Figure 4a). 
According to Table 1, by increasing the angled 
mouth of slots, the number of slots and the 
rotational speed of rotors, the weight of solids, 
would be increased (Figure 4b). In other 
words, the factors chosen in Equations (2) and 
(3) could be verified. Additionally, the effect 
of each treatment and the dual interaction 
between them were significantly higher for the 
weight flow rate of solids (Table 2), and the 
coefficient of estimate of each factor was 
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shown in Table 3. According to Table 3, since the effect of those factors would be examined 
the coefficient of estimate of single effects was by the RSM method simultaneously while the 
greater than dual effects, the results were corn seed specifications were examined by the 
presented based on the single effects. Then, dimensional analysis. 
350 Rl —@— N=4 ——— N=5 N=6 
E 300 _ 300 
: o 
eb R4 ‘a 250 
2 250 a E 
= BR? > 200 
= = 150 
2 150 R3 = 
Ss 
RO 
& RO % 100 
= 100 = 
= s 50 
50 ZB 0 
39 43 47 51 55 59 63 67 71 75 79 22 24 26 28 
n (rpm) Angled mouth of slots (deg) 


(a) (b) 
Fig.4. Weight of solids versus (a) Rotational speed of rotors and (b) Angle mouth of slots 


Table 2- Analysis of variance of the effect of the treatment on weight flow rate of solids 


Treatments df Mean Square F-Value 
Angle mouth of the slot (B) 2 46075.9910 786.41" 
Number of slot (N) 2 65024.2958 1109.81" 

Rotor Speed (n) 3 27279.9565 465.6 

Nx 4 1865.9782 31.85" 

Bxn 6 1084.8282 18.52" 

Nxn 6 1032.1254 17.62" 

BxNxn i2 106.6767 1.82" 

Model 35 158.72 

Total 108 
cv 4.504 


**= P-value < 0.01 , ns = not-significant 


Table 3- Coefficient of estimate of each factor 


Factor Coefficient of estimate 
Angle mouth of the slot (B) 742.48 
Number of slots (N) 865.11 
Rotor Speed (n) 758.26 
NxB 176.26 
Bxn 258.83 
Nxn 246.05 
Be 90.39 
N’ 300.6 
n -230.43 


Intercept 3386.48 
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Development of model, dimensional analysis 
approach 

In the first step, 80 percent of the data was 
randomly selected to develop the desired 
equation. As mentioned, the data collected 
were categorized based on _ dimensional 
analysis according to Equation (2). Initially, 
logarithms of both sides of Equation (2) were 
used to have linear correlation, and then 
equation groups were correlated as previously 
worked by Jafari (1991). 

Log(nt?Q~*) 

= LogC 

+K [Loge (envi) + Log, (en e*) + Log@3(th7*) (6) 
+ Loge, (NtD.*)+Loges(B)+Loge,(Del-*) 

+ +Logp,(Ded-)| 

In Equation (6), C and K represent the 
intercept and the slope of the logarithmic line, 
which can be calculated after founding all 
functional correlations of the above equation. 
Therefore, for convenience, log C and K were 
considered to be O and 1, _ respectively 
(Karparvarfard and Rahmanian-Koushkaki, 
2015). 

Log(nt?Q"*) 

= Logg, (ents) + Log, (ernie) + Loge3(th +) 
+ Loge,(NtD.*)+ Loges(8)+Logep,(Del~") 

+ +Log@,7(D.d~) 

The functional g; through @7 would be 
analyzed in the same order based on the first 
residuals to sixth residuals definitions. For 
example, the first residuals, as an independent 


(7) 


term, represent the values of Log@, (en e) 
(Equation 8) (Jafari, 1991) which has been 
subtracted from the corresponding value of the 
left-hand side, as dependent term, of equation 
(6) to remove the effect of this function (@ 1) 
and prepare to set up the next function. 

Then the first residual (Log@, (ene) 
would be plotted versus their corresponding 
values of exp(n,/tg~+) in this effect, and 
Log@, (ery te) would be obtained (Equation 
9). This order would be repeated for second 
residuals ((log@3(th™*)) that means 
subtracting the values of Log@, (env te) 


from the corresponding values of the first 
residuals in order to remove its effect and 


would be plotted versus their corresponding 
values of the log(th™*), and the log@3(th™*) 
would be obtained (Equation 10). 

Functional relationships for g; through @7, 
have been already obtained by (Karparvarfard 
and Rahmanian-Koushkaki, 2015) and were 
recorded here as below: 


loge, (ev) 

= [19.902 log(th~*) (8) 
+ 27.093 log(NtD.~*)] exp (n/tg-*) 

loge, (e™v'8*) = —6.728 exp(n/tg-) (9) 
log; (th-1) = —14.092 log(th-*) (10) 
logp,(NtD, +) = 5.2324 log(NtD, *) (11) 
log ps5(B) = 8.6585 log(B) (12) 
logs (Del!) = —34.996 log(D.1"*) (13) 
log; (Ded~*) = 34.752 log(D.d“*) (14) 


At this stage, equations (8) to (14) could be 
combined with Equation (2) as follow: 
log(nt?Q-1) = logc 

i K{[19.901 log(th~?) 

+ 27.093 log(NtD~) 

— 6.728]envie* 

— 14.092 log(th~+) (15) 
+ 5.2324 log(NtD, *) 

+ 8.6586 log(8) 

— 34.996 log(D, 1-1) 

+ 34.752 log(D.d-*)} 


The final value of K and log C were 
obtained by plotting the corresponding value 
of log(nt?Q~*) against the numerical values 
of the expressions inside bracket of equations 
(14) (Figure 4). 

K = —0.0352 and logC = 0.6593 

The final equation was concluded: 

nt?Q°? 

= 4.564 {[cr-p recon, 
= 14.092 

10-79] E(B)  ea,y92" (16) 


x (8) 6845 x (D171) 34-973 
—0.0352 


= ae 


(1.87 


x (D a 
e 


It should be noted that the present study 
was carried out within the limits of the values 
in Table 4. Therefore above equation can be 
used for predicting weight flow rate discharge 
of a multi-slots feeding device of a row-crop 
planter within the limits mentioned in Table 4. 
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0.2 
y = -0.0352x + 0.6593 


R?= 0.9491 


log(nt’ Q ~') 


12 14 16 18 20 22 
log (p_1*@_2*p_3*p_4*@_5*p_6*Q_7) 


Fig. 5. Plot of log(nt?Q*(—1)) versus log(@1- Pz. 03. Pa: Ps. Yo P7) 


Table 4- Considered a range of dimensionless groups in the present study 


Variable groups Range of variations 
th-? 1.646 — 2.073 
NtD,* 1.226 — 1.913 
n/tg-? 0.0269 — 0.0525 
B 23 <27 
DI? 5.004 — 5.206 
D.d7t 6.460 — 6.721 
Development of model, response surface RSM modeling. The result was the following 
methodology approach equation which can be used to predict the bulk 
The purposing of modeling performance for volume flow rate of solids. 
the single-seed metering device is crucial to Q = 31790.1 — 1604.92 x B — 5115.55 
establish the correlations between the XN — 7245.58 xn 
requested inputs (for instance: angle mouth of + 88.1321 xB xN 
slots, number of slots and the rotational speed eal XB xn 
+ 787.949 xNxXn 
of rotor) and outputs (bulk volume flow rate of + 22.5983 x B? + 300.6 
solids). The RSM as a statistical model was +x N2 — 2363.13 xn? 
applied for each performance parameter of the R’ = 0.96 


planting device using Design-Expert software The ANOVA table of the RSM model is 


(Version 11.0.1). Table 5 shows the maximum presented in Table 6. 
and minimum of data which were used in the 


Table 5- The maximum and minimum of data sets used for RSM modeling 


Variable Maximum Minimum 
B (Degree) 25 23 
N 6 4 
n (rps) 1.29 0.67 
Q (mm?.s") 7025.029 1639.235 
Table 6- Anova Table of RSM model 
Source Sum of square df Mean square __F-Value 
Model 3.558E+6 3 1.186E+6 26.57" 


**= P-value < 0.01 
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Accuracy analysis of the RSM and dimensional 
analysis models 

In order to assess the validity and accuracy 
of estimating the volume flow rate of seeds by 
the two proposed models, 20 percent of the 
data was set aside to be used for verifying the 
two models. The volume flow rate as predicted 
by dimensional analysis and the RSM model 
were plotted against the measured volume 


® Dimensional Analysis 


flow rate (Figure 6). The plot shows that the 
values predicted by the dimensional analysis 
are closer to the measured values as compared 
to values predicted by with the RSM values. 

The standard error of estimate was 68.13 
mm? s’ and 475.59 mm* s’ for the 
dimensional analysis and the RSM method, 
respectively. 


RSM 


_~ 7000 
Yr Q (Dimensional Analysis model) = 1.0052x - 29.378 
on 2— 
: 6000 R2 = 0.9963 
~~ Q (RSM model) = 0.7828x + 633.07 
g 5000 R2 = 0.864 
heal 
3 
lo) 
- 4000 
= 
a 
© 3000 
> 
lao} , 
i) f 
2 2000 
mo) 
(0) 
a 
® 1000 
& 
onl 

0 

0 1000 2000 


3000 


4000 5000 6000 7000 


The measured volume flowrate (mm? s”!) 


Fig.6. Plot of predicted volume flow rate versus measured volume flow rate 


Conclusions 


In this study, two equations were developed 
which can be used to predict the volume flow 
rate discharge of corn seeds from a new type 
of planter metering device the following 
conclusions were also made: 

1. The increase in rate for a bulk volume 
flow rate of solids showed a descending slope 
by increasing the rotational speed of rotors. 

2. More effective performance for bulk 
volume flow rate of solids can be achieved at 
lower rotor speeds. 
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